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I. INTRODUCTION 



The classical baker map on the torus is often used to illustrate properties of chaotic dynamical systems . It 
may be regarded as the simplest generalization of the Bernoulli shift in two dimensions. The baker map is measure 
preserving and chaotic: it belongs to a small class of the dynamical systems, for which the metric Kolmogorov- Sinai 
entropy is positive and may be calculated analytically. This classical dynamical system has been generalized in several 
ways. The generalized baker map was used to study properties of invariant measures and fractal dimensions 
while snapshot attractors appearing in the random baker map were analyzed in . Another generalization of the 
model, called multibaker map |^,^ , is useful for investigating the effects of deterministic diffusion. Chaotic dynamics 
of a classical map can be described as a linear evolution of densities transformed by an associated Frobenius-Perron 
operator. Much progress has recently been achieved ||8|,[l0|-|l^ in finding the spectrum and understanding some unusual 
properties of the right and the left eigenstates of the F-P operator of the classical baker map. 

A first attempt to quantize the baker map on the torus is due to Balazs and Voros fl^ ]. Another, more symmetric 
version of the quantum map was found by Saraceno . A general construction allowing one to quantize a class of 
piecewise linear maps on the torus (including the baker map) was given by De Bievre, Degli Esposti, and Giachetti 
While another quantization scheme was recently proposed by Rubin and Salwen p8[ | . The quantum baker map on the 
torus become a standard model used to pursue the concept of semiclassical quantization of non-integrable classical 
systems p^- p5| , p7[ -p9| emerging from the theory of Gutzwillcr |3^j3l[] . In several other papers devoted to the quantum 
baker map one demonstrated sensitivity of ei gen states of the system with respect to small perturbations [ ^2[ , proposed 
an optical realization of this quantum map ||33|1 , and designed a scheme of realization of the baker map by methods 



of quantum computing 34 



In this work we propose a different dynamical system: the baker map on the sphere. In fact we define an entire 
family of classical systems and the corresponding quantum maps. On one hand, the classical version of this model 
is, to our best knowledge, the first dynamical system on the sphere with explicitly computable, positive K-S entropy. 
On the other, analysis of several possible versions of the analogous quantum systems may help in clarifying various 
aspects of the quantum-classical correspondence for non-integrable dynamical systems. 

In the simplest case of the model the map defined on the sphere is not symmetric with respect to the reversal of 
time. The Floquet operator describing the time evolution of the corresponding quantum map is constructed out of 
the Wigner rotation matrices and is represented by orthogonal matrices of even dimension. In the general case the 
evolution operators are represented by complex unitary matrices. 

Each quantum model, corresponding to a given classical system, can be generalized by introducing into the model an 
arbitrary complex phase [ p^ . One obtains than a one parameter family of quantum maps acting in the same Hilbert 
space, which may be useful for studying various aspects of level dynamics: curvatures, correlations of velocities and 
anticrossings. 

Phase space features and the time evolution of the quantum map can be suitably analyzed with help of the SU (2) 
coherent states. A simple generalization of the quantization scheme presented in provides four different maps 
(with different spectra and traces), that seem to correspond to the same classical system. This quantization scheme 
allows us to obtain two versions of the symmetric quantum baker map on the torus different from that introduced 
inll]. 

The paper is organized as follows. The quantum version of the baker map on the sphere is introduced in Section 
II. The corresponding classical system is analyzed in Section III, while the statistical properties of the quantum map 
are discussed in Section IV. Appendix |^ contains the proof of unitarity of auxiliary matrices constructed out of the 
Wigner rotation matrices. Detailed formulae for four versions of the quantum baker map on the sphere are given 
in Appendix while the entire family of the classical and quantum models is given in Appendix |c[ An analogous 
construction for the baker map on the torus is shown in Appendix 

II. CONSTRUCTION OF THE QUANTUM MAP 

To construct the quantum map corresponding to a given classical transformation on the sphere we use the iV- 
dimensional representation of the rotations group. Using the angular momentum operators commutation relation 

[Ji, Jj] — ieijkJk (1) 

we obtain N = 2j + 1 base vectors \j, m)i, eigenstates of and Ji (i-th coordinate of J) 

J^\j,fn)i ^ j{j + l)\j,m)i, J,\j,m), ^ m\j,m)i, (2) 
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m = —j, ■ ■ ■ ,j. To work in an even dimensional Hilbert space Ti. we consider the half-integer values of j. 

We design the quantum baker map Bs on the sphere in such a way that Bs commutes with and the quantum 
number j is preserved. Therefore it is convenient to simplify the notation by writing \m)i instead of \j,m)i. The 
action of the baker map has to stretch phase space along one direction (z-axis) and to squeeze along the other one 
(say, along the x-axis). The transformation between these two basis \'^)x — R\'^)z is described by the Wigner rotation 
matrix R representing a rotation around y-axes by the angle tt/2. We use its representation in the |m)z basis 

R,n',m^ z{m\e-'^^y\m'), (3) 

with m, m' — —j, This matrix is real and has the following property 

p — 1"!™^™ ?? ( 1 F? (A\ 

Let us split the Hilbert space TL into two subspaces distinguished by its projection on \m)z 

n^' = (l^-) e n : ,{m\'^) =0, m > 0} 

= (l^-) g H : ^(ml*) = 0, m < 0}, (5) 

so that contains the states localized on the southern hemisphere and - on the northern one. Projection on 
Ifc)^: leads to another partition 

H^ = {|^r) =0, fc>0} 

7^^-{|*> e?i :,(fc|*> -0, fc<0}, (6) 

soT-L = ^ T-C^' + n^. The baker map transforms each vector j^''^') e Ti'^ into a vector 1$^) e Ti^ , 

with a linear stretching by a factor of two in the z direction. The effect of this stretching can be formalized by the 
condition that the N/2 odd (even) coefhcients of the state |<i>^) expanded in the basis \m)z are given by the first 
N/2 coefficients of the state l^*"^) represented in the same basis 

,{m\-^^) ^ V2 z{2m + for to = -j, . . . , -1/2. (7) 

The factor V2 leads to the correct normalization. 

The analogue condition for the linear squeezing in a;-direction would give the overcomplete system of equations with 
no solution. This fact is rather intuitive because it is not possible to map the southern hemisphere into the western 
one with the linear stretching in z-direction and the linear squeezing in x. However, we require |<i>^) to be squeezed 
to the western hemisphere H,^ 

,(A;|$^)=0 for fc-l/2,...,j. (8) 

The above condition corresponds to a nonlinear squeezing in a:. We want to find a linear operator B^^ which 
transforms into |$^) 

\^^)^ = B''^\^').. (9) 

To rewrite this formula into the \k)x basis we multiply its both sides by the transformation matrix R: \^^)x = 
RB^^ \'^^) z- To satisfy the condition (||) all elements of the lower half of the matrix RB^^ must vanish. Since we 
want RB^^ to act only on the subspace Ti.^ we set to zero the right half of this matrix. Denoting by M the only 
unknown N/2 x N/2 block of the matrix RB^^ we put equation (^ in the form 



' ' Afi,m 



l=-j m=-j 







toI^-^) , (10) 



Inserting above into (|^) we obtain the following conditions for the matrix M 

_ j_ 

^E 2fc+,;^^'^™ = '^fc,™ for fc,TO = -j,...,-i . (11) 



In the appendix ^ we show that the matrix R' composed of half of every even (odd) column of the matrix \/2 R is 
unitary, (since it is real, it is also orthogonal). Consequently, we take Mi^m — Ri,2m+j ='■ R'l m (': = ■ • 1/2) 
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The other matrix i?"^ 



N 



and insert it into ( p^ constructing in this way the transformation : Ti.^ — > Ti.^ 

Ri,2m-j (now Z, m = 1/2 . . . j) is also unitary, which allows us to build the second transformation B^^ : TL 
in the same manner. We define the quantum baker map on the sphere as a linear operator Bs = B^^ + B^'^ 
on 7Y. In the \m) ^ basis the quantum baker map on the sphere takes the form 



R' 
R" 



acting 



(12) 



In the mixed representation {z{'m\Bs\m')x) the baker map on the sphere consists of two diagonal blocks R' and 
R" . Figure |l| shows (Bs)^ in this representation for n — 0,1, 2, 3. The first plot obtained with n = displays just 
the elements of the rotation matrix Rm'm = z{n^\TT^')x- For m and m' fulfilling + m'^ > such elements are 
exponentially small. This fact has a simple classical analogy. Consider a partition of the unit ball into N slices of the 
same width perpendicular to the direction z and another partition consisting of N slices perpendicular to the direction 
X. Some side z-slices (of small radii) do not overlap with some side x-slices. This 'circular' structure manifests itself 
also in the pictures obtained for n = 1, 2 and 3 and contrasts with the rectangular patterns visible in the analogous 
figure drawn by Saraceno and Voros pSl for the quantum map on the torus. 




FIG. 1. Quantized baker map on the sphere (-Bs)" shown in the mixed representation for TV — 128 and n = 0, 1, 2 and 3. 
Gray tone scale represents values of real elements of TV x TV matrices. 



The rotation matrix R is unitary, as are the matrices R' and R" , constructed out of its elements. Hence the quantum 
map -65 is unitary. Writing equations (|^) describing the double stretching of the vector we can choose ^(2m + j| 
(odd column of R) or z{2m + j + 1\ (even column) projections of the state |$^). A similar ambiguity remains in the 
choice of the vector |$^), so one may construct four different versions of the quantum map: Sg*™', Bs*°^\ Bs*^"' 
and -65^^^^ as written explicitly in Appendix In the classical limit TV —s- 00 all four families of quantum maps seem 
to tend to the same classical system. 

In order to create one parameter families of quantum maps it is possible to introduce a phase factor A into any of 
the above four versions of the model 

■'^R' 
e~'^R" 



BsiX) = R- 



(13) 



This particular way of introducing the parameter A is advantageous, since there is no drift of the eigenphases (i.e. 
the phases of complex, unimodular eigenvalues of Bs) with the parameter (the mean velocity {d(j)i/dX), averaged over 



4 



individual eigenphases, vanish). An analogous generalization of the quantum baker map on the torus was already 
proposed by Balazs and Voros [p^ . It gives a family of quantum maps corresponding to the same classical system. 

Moreover, it is possible to generalize the model by changing the direction of squeezing, which corresponds to a 
different choice of the rotation matrix R. As discussed in Appendix |c[ these models lead to a family of the classical 
systems, which can be parametrized by the angle 7. The classical system corresponding to 7 = 7r/2 (rotation around 
the X axis) possesses a generalized time reversal symmetry, so varying the parameter 7 one can study the effect of the 
time reversal symmetry breaking. 



III. CORRESPONDING CLASSICAL MAP 



The quantum baker map on the sphere is unitary, hence the corresponding classical map M : S"^ ^ S"^ has to 
conserve the volume of the phase space. The quantum transformation corresponds to the stretching by a factor of 
two along the z axis. Thus it is linear in the variable t — cosO. The operation of squeezing is nonlinear in a;, but must 
be linear in the angle which leads to the following classical map on the sphere 

{{21 - 1, ip/2) ioT t > and if < TT 

{2t + l,ip/2 + TT) for < < and < TT , . 

{2t - 1, ip/2 + tt) for t > and ^ > TT ' ^ ' 

{2t + 1, ip/2) for i < and <^ > TT 

It is presented in Fig. |^ in the (t, ip) coordinates. 




FIG. 2. Classical baker map on the sphere in the coordinates {t = cos6,ip). Parts Ai, A2, B\, B2 are linearly transformed 
into A'l, A'2,B[,B'2. 

The classical baker map on the sphere is chaotic: the dynamical entropy of Kolmogorov-Sinai equals to In 2. The 
generating partition consists of two cells A — Ai (J A2 and B — Bi U B2 as shown on Fig. ||. The uniform measure is 
invariant under the map ( p^ and after each iteration all probabilities Paa, Pab, Pba, Pre of going from one cell to 
the other are equal. A similar reasoning done for n time steps shows that the probabilities of all possible 2" trajectories 
are equal to 2~", which leads to the same metric entropy as for the baker map on the torus. A generalization of the 
model leading to a family of the classical systems characterized by the same dynamical entropy, is given in Appendix |^. 

Furthermore, the topological entropy (as all generalized Renyi entropies [p6|) is equal to In 2. Therefore, the number 
of points belonging periodic orbits growth with their period n as 2". In particular the number L of periodic orbits of 
(BsY^ is equal to 2" for n odd and 2" — 2 for n even. In analogy to the work of Saraceno and Voros |2^ we constructed 
a generating function which allows us to find the action of each trajectory, required for the semiclassical treatment of 
our model |37[] . 

In order to demonstrate a correspondence between the classical and the quantum baker maps on the sphere we use 
the vector coherent states ^8|-^. Each point on the sphere, labeled by the spherical coordinates {0,ip), corresponds 
to the SU{2) coherent state \9,ip), defined as 



1^,1^9) = exp i9(s'mipJx — cos (pJy) 

Expectation values of the components of the angular momentum operator J are 

{9, ip\J\9, if) — j(sin9cosip, sin0sin(/3, cos^). 



(15) 



(16) 



which establishes the link between the coherent state \9, ip) and the vector {9, ip) oriented along the direction defined 
by a point on the sphere. Simple expansion of vector coherent states in the \m)z basis |^ , ^ 
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exp — 'm)ip) 



2j 



1/2 



(17) 



makes them handy to use in analytical and numerical investigations of any quantum system corresponding to a 
classical map on the sphere. 

Correspondence between the classical and the quantum model is pointed out in Fig. ^. The right column presents 
contours of the function 



FU0,cp) = \{e,^\{Bsr\9,cp)\' 



(18) 



labeled by the number of iterations n. Peaks of this function correspond to the states for which the probability of 
staying after n iteration of quantum map is maximal. The left column shows periodic points of the classical map 
(p^). The north and south poles (represented in Fig. ^ by horizontal lines) are fixed points of the classical map, but 
the map is not continuous there. All pictures are plotted in the (t, if) coordinates. Quantum data are obtained for 
N — 200 and the variant -Bs'"^-* of the model (see appendix]^). The analogous pictures done for other variants of the 
quantum map look the same. 




FIG. 3. Comparison between the classical periodical orbits of length n (crosses in the left column) and squared amplitude 
of quantum propagator \{6,if\{Bs)"\d,tp)'^ (right column) for the baker map on the sphere. There exist L periodic points of 
(Bs)" (crosses). Two of them are localized on the north and south poles, and are symbolically represented by stars arbitrarily 
drawn at (p = n. 

Analysis of eigenvectors \vi) of the quantum map Bs provides an additional support for the quantum-classical 
correspondence. Generalized Husimi distribution of some eigenvectors, Hv-(9,(p) = |(t;i|6', (p)\'^, reveals maxima in the 
vicinity of some classical periodic orbits. This effect, called quantum scars Gsl, was already observed for the standard 
baker map on the torus [[l6| . 

The correspondence between a given classical system AI : X X and a family of quantum system Uj : Ti. ^ Tl, 
parametrized by the quantum number j, may be quantitatively characterized by the following condition of regular 
quantization. Let C[{6, ip), p] denote the circle on the sphere of radius (angle) p centered at {d,ip)- Let us define the 
quantity 
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i^{e,^):^^i-tl f \{e',^'\u,\e,^)\^sine'd9'd^' (19) 

JC[M{e,ip),p\ 

measuring localization of the Husimi function of the transformed state Uj\9, ip) in the p-neighborhood of the classical 
image M{6, Lp). Due to the normalization factor (2j + 1)/Att the integral of the Husimi function over the entire sphere 
X is equal to one. If for any p > 

lim inf [Ipie.'f)] ^ 1, (20) 

the quantization procedure linkingthe classical map M and the family of quantum maps Uj is called regular with 
respect to SU{2) coherent states Q,|5|. The infimum is taken over all points {9, ip) of the classical phase space, and 
the size of the Hilbert space A'^ = 2 j + 1 serves as a parameter in the family of quantum maps Uj . 

Not being able to prove this condition analytically, we took for M and U the classical ( |T^ ) and the quantum ( |T^ ) 
baker maps on the sphere, respectively, and performed extensive numerical tests studying the dependence of Ip on N 
and p. Figure |4| shows how the quantity Ip tends to unity in the semiclassical regime of large N. The mean value of 
Ip, averaged over 100 points on the sphere, is close to unity at A'^ = 600 (o), but the minimal value converges much 
slower (x). In spite of this fact, the above results confirm the existence of a tight relation between the classical and 
the quantum models of the baker map on the sphere, introduced in this paper. 








200 400 

N 



600 



FIG. 4. Regular quantization condition linking the classical ( [L4|) and the quantum ( |l2| ) baker maps on the sphere. Depen- 
dence of the locally integrated Husimi functions Ip on the size of the matrix A*' for p — n/lb. Symbol (o) represents the mean 
value averaged over 100 points placed uniformly on the sphere, while the maximal and the minimal values are denoted by (+) 
and (x), respectively. 



IV. STATISTICAL PROPERTIES OF THE QUANTUM MAP 



Quantum maps corresponding to classically chaotic systems are expected to display the spectral fluctuations char- 
acteristic of random matrices pq , ^ . Due to the lack of the time reversal symmetry for the baker map on the sphere, 
unitary matrices Bs should be compared to the Dyson's circular unitary ensemble (CUE) E^,^. To obtain a sat- 
isfactory statistics we accumulated data from all four versions of the model (see appendix^) and several sizes of 
matrices. 

Eigenphases (f>i of random unitary matrices are distributed uniformly in [0,27r], so no unfolding of the spectrum is 
necessary. We numerically diagonalized matrices Bs, ordered the eigenphases and computed the rescaled spacings 
si — (0/+1 — (t>i)/ (s), where / = 1, . . . , iV — 1 and the mean spacing (s) is equal to 2'k/N . Figure || presents the near 
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neighbors distribution P{s) calculated from matrices Bs with dimensions ranging from TV = 50 to iV = 800. The 
dashed line is the Wigner distribution P„(s) = 327r~^s^ exp(— 4s^/7r^) (exact for 2x2 Hermitian matrices of the 
Gaussian unitary ensemble), which provides a very good approximation for the asymptotic CUE result obtained for 
N ^ oo . Collecting more data out of all four variants of the map we observed that the precision of the Wigner 
surmise is not satisfactory any more. Numerical data conform well to the exact CUE result, implemented as a power 
series [Q, and shown in the inset. Note that for the unitary ensemble the differences between the Wigner surmise 
Puis) and Pcue{s) do not exceed 5 parts per thousand. 




FIG. 5. Eigenphases spacing statistics for quantum baker map on the sphere. Histogram contains 79524 data obtained 
for A'' = 50, . . . ,800. Dashed line represents the Wigner distribution Pu{s). Histogram in the inset consists of 318096 data 
and exhibits the difference from the Wigner surmise P{s) — Pu{s). The solid line represents exact results for random matrices 

PcUEis) ~ P^S). 

To analyze the long range correlations we computed the spectral rigidity A3 , introduced by Mehta and Dyson . 
Numerical data obtained from matrices with dimensions = 400 . . . 500 are compared on Fig. || with CUE results 
represented by a dashed line. In addition we verified that the statistics of eigenvectors of matrices Bs fit for large N 
to the x^2 distribution predicted for CUE 

Since the matrix Bs is orthogonal and its spectrum consists of two replicas of the same sequence, only the first N/2 
spacings out of each matrix diagonalized were used in the statistics. Although the joint probability of eigenphases for 
random orthogonal matrices differs from this characteristic of CUE |^^, the level spacing distribution are the same 
in the limit of large matrices. This fact, observed numerically |53| , can be proved in a rigorous way [Q. Therefore it 
is not surprising that the orthogonal matrices representing the quantum baker map on the sphere exhibit CUE - like 
spectra. 

Using SU{2) coherent states it is possible to visualize the propagation of a wave packet. As before we use Husimi 
phase space representation. Fig. ^ presents contours of Husimi functions with <i> = {Bs)"\Oq, (po) drawn in the 
[t, if) coordinates. In other words each graph represents the subsequent iterate of the coherent state initially localized 
at {6o,ipo). After each step the wave packet is squeezed along the (^-direction and stretched along the t-direction 
as predict the formula (14). At n = 4 the state occupies the northern and the southern hemispheres so it is split 
into two part after the next iteration. Dimension of the Hilbert space N is 200. Observe an abrupt change in the 
shape of the wave packet occurring after 7 — 8 iterations. This number corresponds well to the logarithmic time scale 
n* = log2 N 7.6, which determines the behavior of the quantum baker map on the torus [M and often emerges in 



several problems of quantum chaos (see e.g. 
{N = 10,24,50,100). 



55[). The same time scale n* was observed for all other cases investigated 
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FIG. 6. Spectral rigidity for the quantum baker map on the sphere. Crosses denote numerical data while the line represents 
the theoretical prediction for CUE. 

Let us also mention that a generalized version of the model, discussed in Appendix covers also a system enjoying 
the time reversal symmetry. In this case, obtained for the parameter 7 set to 7r/2, the spectral statistics of the Floquet 
operators coincide with the predictions of the circular orthogonal ensemble (COE). 



In this paper we introduced classical and quantum versions of the baker map on the sphere. While the classical 
dynamics takes place on the sphere indeed, the name quantum baker map on the sphere cannot be treated literally: 
the quantum dynamics takes place in a finite dimensional Hilbert space 7i, and a link with the classical dynamics on 
the sphere can be achieved with help of the SU{2) coherent states. Nevertheless, our quantum model differs, in many 
respects, from the quantized baker maps of Balazs and Voros [|l5| and Saraceno In a specific case of the model it 
is represented by orthogonal matrices, which exhibit CUE ~ like statistical properties of the spectra. This fact reflects 
the lack of any antiunitary symmetry in this system. Since a time reversal case of the model, displaying a COE - like 
spectra, was found, the model proposed may be used to study the effects of the breaking of a generalized antiunitary 
symmetry. 

Our construction is based on the properties of the Wigner matrix, representing rotation by the angle tt/2. For each 
value of the classical parameter 7 we found four different families of quantum systems, parametrized by the even size 
iV the Hilbert space, which in the limit N ^ 00 seem to correspond to the same classical system. The condition of 
regular quantization with respect of SU{2) coherent states has been checked to a satisfactory precision. 

In the simplest case of the model (7 = 0), the quantum baker map on the sphere is represented by real matrices. Thus 
the imaginary part of their traces vanish, in contrast to the quantum maps on the torus which suffer a logarithmical 
divergence of the imaginary part of the trace in the semiclassical limit |^ . 

If Fourier matrices are used in our construction in place of the Wigner rotation matrices, we obtain four versions 
of the quantized baker map on the torus, some of them different than these previously known. We hope that the new 
dynamical systems will prove their usefulness in further studies on the quantum - classical correspondence for chaotic 
systems. 



V. CONCLUDING REMARKS 
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FIG. 7. Evolution of the coherent state localized at {^tt, ^tt) under action of the baker map Bs for A'^ = 200. Contours 
of Husimi function are plotted in the {t, ip) coordinates. Note the delocalization of the wave function for n > log2 200 ~ 7.6 
iteration. 
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APPENDIX A: ROTATION MATRIX R AND UNITARITY OF R' 



Components of the N dimensional Wigner rotation matrix R, describing rotation around y-axes by the angle 7r/2, 
can be written as [351 



= 2- V(J + -m) ! {] - m) ! (j + k) ! (j - k) ! Su,m , 
where N — 2] ^ 1 \s even and k^m — —7, —j H- 1, . . . , j — 1, j. The matrix Sk,m reads 

min(j+fc J— m) 



k.7n 



E 



p— max(0,A:— m) 

We want to show the following identity 

R-k,m = {—^y~™'Rk, 



p\ {p + m ~ k) \ {j — m ~ p) \ {j + k — p)! 



(Al) 



(A2) 



(A3) 
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To this end we start computing R-k,; 



R^k,m = 2-V(j + mV-U - my.{j ~ ky.{j + ky. s-k,m ■ (A4) 

Rearranging the sum we get 

^ p\{p + m + ky.{j-m-py.{j-k-py.' ^^^"^ 

p— max(0,— fc— m) 

In ( [A5| ) we replace phy q = j — m~ p (so p = j — m — q) and determine the new hmits of summation: qmin = j — m — 
Pmax — —TTi — min(— fc, — m) = max(0, A; — m) , while (/max = j Pmin = J — m + min(0, fc + rn) = min(7 + fc, j — rri) . 
Now we rewrite the sum (A5) 

" 5* ^ (j-m-g)!(j + fc-g)!g!(g + m-fc)! ' ^^^^ 

g— max(0,A: — m) 



Comparing the result with (JA^) we can notice that 

S-k,m - (-l)^-'"5fc,m (A7) 



and the identity (A3) follows. 

Unitarity of the N/2 x N /2 matrix R' will be deduced from the unitarity of R 



E ^Vfc^fe,™ = E ^fc.'^fc,™ = -5;,™ ■ (A8) 



This sum may be divided into two parts 



E {^kdRk,m + R-k,lR-k,m) = Sl^m ■ (A9) 
k=-j 



Now we use property (A3) to reformulate the first part of the sum 

_ i_ 

E Rk,lRk,rn (1 + (-1)2^-'-'") = 6l,^ . (AlO) 

Choosing odd columns of matrix R reduces to exchange I 21 + j, and m 2to + j, so 2j — / — m — > — 2(/ + to). 
Since N is even, j and the indices I and to are half integer, so Z + m is an integer number. We can rewrite the result 

2 

2 E Rk,2l+jRk,2m+j = ^2l+j:2m+j = ^l,m ■ (AH) 
k=-J 

This proves that choosing upper halves of even columns of the matrix \/2R gives an unitary N/2 x N/2 matrix R' . 
Choosing odd columns or bottom halves {k > 0) gives other unitary matrices of the size N/2. 



APPENDIX B: FOUR VERSIONS OF QUANTUM BAKER MAP ON THE SPHERE 

As shown in appendix ^ the auxiliary matrices R' and R" are unitary, and can be constructed out of the elements 
of the rotation matrix R in several ways. Let us introduce the indices a, 6 = 0, 1 denoting, whether the odd or the 
even columns of R where used in the construction. This allows us to define four unitary matrices 

:= V2Rk,2i+j+a, fc, / = -J . . . - i; R"f} V2Rk,2i-j-i+b, kj^^...j; (Bl) 
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which lead to four different variants of the quantum baker map on the sphere for each even N 



Bs 



(ab) ^ ^-1 







(B2) 



with a, b — 0,1. Each of matrices Bs^""^^ is orthogonal and my be generalized into a one parameter family of unitary 
matrices according to formula ([T^). Fo r co ncreteness we provide here some examples. For N = 2 any unitary matrix 
is equivalent to rotation, and formula (BS) gives 



V2 



1 1 

1 -1 



Simplest interesting matrices appear for iV = 4 



and Bs 



(01) 



1 -1 



1 

V8 



2 





2 





V3 


1 


-V3 


-1 





2 





2 


-1 


V3 


1 


-V3 





2 





-V3 


1 


1 


V3 


1 


2 





V8 





2 


-1 






_ -1 


x/3 





2 



1 



x/3 


-1 


2 





2 







-1 


1 


V3 





2 





2 


1 


-V3 



B: 



,(11) 



1 

7i 



rv3 


-1 


-x/3 


1 


2 





2 





1 


V3 


-1 







2 





2 



Observe that the first N/2 columns of Bs(°°) and Bs'-°^^ are equal, while the same is true for Bs*"' and Bs'-^"\ 
Moreover, the matrices Bs*°°' and Bs''^°^ share the same last N/2 columns; the same is true for Bs''^^'' and Bs'^^-*. 
These relations, valid for arbitrary even N, follow directly from the definition (B2). To emphasize the differences 
between all four variants of the quantum system let us consider the trace of each matrix. It is easy to show that for 
any N the matrices ^5^°' are traceless: TrBs'^°' = 0. Moreover, TiBs^°^^ = V2, while TrBs''°°^ +TrBs^^^^ = V2. 
For = 6 the traceless baker map is represented by the matrix 



RjlO) 



8V2 



3V5 -VTO 3 
8 

3V2 6 -VTo 

8 

-1 3V2 3%/5 
8 



8 
3V5 -3V2 

8 
/To -6 


-3 VTo 




1 


-3%/2 
8 

-3\/5 



APPENDIX C: A FAMILY OF SYSTEMS PARAMETRIZED BY A CLASSICAL PARAMETER 7 

Discussed quantum system may be generalized by picking for R a different rotation matrix. Let us allow for a 
rotation along an arbitrary axis, which belongs to the the plane xy and forms the angle 7 with the y-axis. The 
transformation matrix from the \m)z basis to the new basis takes the form 



(CI) 



Replacing the rotation matrix R by in the formulae (Bl) and (B2) we obtain a continuous family of quantum 



models B.^"''' parameterized by the angle 7 € [0, 27r) and corresponding to a various directions of squeezing. 
We can find the classical counterpart of the generalized quantum model 



2t- 1, 
2t+l, 



TT 
TT ■ 



27) mod 27r] /2 + |7r + 7) mod 27r^ for t > 
27) mod 27r] /2 + Itt + 7} mod 27r) for t < 



(C2) 



This reduces to (^4|) for 7 = 0. Figure || shows how the sphere is mapped onto itself by the formula (|C^). The sphere 
is cut along meridian (</? = tt + 27). The northern and the southern hemispheres are stretched in the t-direction by 
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the factor of 2 and are squeezed in the tp-direction. The transformed blocks are respectively placed on the east and 
on the west from the meridian {ip = 7r/2 + 7). 




FIG. 8. Generalized classical baker map on the sphere drawn in the {t,ifi) coordinates for 7 = n/6. The parts Ai, A2, Bi, 
B2 are linearly transformed into A'^, A'2, B[, B'2. 

For 7 = 7r/2 the matrix describes the rotation by tt/2 around the x-axis. In this case the classical map is 
symmetric with respect to the reflection along the a;-axis and has a generalized time reversal symmetry. Existence 
of the geometric symmetry corresponds in the quantum case to the fact that the Floquet operator B^yj* commutes 
with the reflection matrix Vnn' = Sn,N-i-n' ■ Therefore one can divide its eigenstates into two classes: the symmetric 
and the antisymmetric states. Analyzing the spectral statistics separately in both parity classes we find the COE-like 
behavior of the spectrum. By varying the classical angle 7 in the vicinity of 7r/2 one may study the influence of the 
time reversal symmetry and the geometric symmetry on the system. 



APPENDIX D: FOUR VERSIONS OF QUANTUM BAKER MAP ON THE TORUS 

In this appendix we obtain four variants of the baker map on the torus and show their relation to the earlier models 
1^,^ . Consider the generalized Fourier matrix p5| which is the transformation matrix from position q to momentum 
p basis on the torus 



1 



-27Ti{k+Xg}(l+Xp)/N 



N 



(Dl) 



The phases XqiXp ^ [0: 1) may be treated as free parameters, which correspond to the phases gained by the wave 
function after translation by one period in the q and in the p direction, respectively. The matrix F^^''^^ is unitary. 
Assuming that the matrix size N is even, we may apply the procedure of taking every second half of column to obtain 
four auxiliary matrices of size iV/2 



■I XqXp (a) 



N/2 



k,l 



f; 



It XqXp (b) 



N/2 



k,l 



k+N/2,2l+b 



(D2) 



where a, 6 = 0, 1, and the indices run fc, Z = 0, . . . , N/2 — 1. Both matrices F' are unitary, since 

N/2-1 N/2-1 

E±_ g-27ri(fc+x,)(2;+a+Xp)/A' g27rj(fe'+x,)(2i+a+Xp)/A' ^ ST^ _ 2TTi(k' -k)(2l+a+Xp) / N 
N ^ N 

1=0 



1=0 



^2TTi{k'-k){a+Xp)/N 



N/2-1 



y J_ e^-^ik'~k)i/(N/2) 
^ N/2 ' 



1=0 



and the same holds for both matrices F" '■''^ . Hence the four versions of quantum baker map on the torus take the 
familiar form 



pl XqXp (a) 
^ N/2 





F] 





/' XqXp (b) 
N/2 



(D3) 



The matrix Bt^^ ^^^^ is equivalent to the original quantum baker map of Balazs and Voros 1 15|, while the symmetric 
map of Saraceno llq] is defined as follows 



13 



B 



Fn 



1/2,1/2 



F 1/2.1/2 





W 1/2.1/2 



(D4) 



Two versions of our map Bt^^^'^^^^^^\ 5^1/2.1/2(01) ^^sq symmetric in respect to reflection {VBtV ^ — Bt 
where Vn,n' = Sn,N-i-n'), as the map of Saraceno. The real part of their traces do not change with N and is equal to 
and ^/2, respectively. They correspond to the same classical baker map on the torus and might be useful for further 
semiclassical investigation of this model. 

Statistical analysis of spectra of these two variants Bt^^^'^/^ 5^,1/2,1/2(10) geveral dimension N e [50,600] 
shows that the level spacing distribution P{s) conforms to predictions of circular orthogonal ensemble with a precision 
allowing one to discriminate the Wigner surmise. The same was checked for the matrix of Saraceno. Due to symmetry 
of the problem the statistical data were collected separately in each parity class containing N/2 eigenvalues. 
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